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The  nonlinear  thermoviscoelastic  constitutive  behavior  of  a  conventional  tank  gun  propellant,  M30,  is  modeled  using  a 
moditied  superposition  integral  that  incorporates  the  effects  of  microstructural  fracture  damage.  This  work  represents  a 
thermoviscoelastic  extension  of  previous  work  (Gazonas  1992, 1993)  that  modeled  the  loading  rate  behavior  of  the 
propellant  over  nearly  five  decades  in  strain  rate.  In  this  w<»k,  uniaxial  compression  tests  are  conducted  on  the  propellant 
at  constant  strain  rates  of  .01  and  210  I/s  and  temperatures  of  -30  and  60  degrees  Celsius.  A  master  relaxation  modulus  is 
determined  from  the  data  in  the  form  of  a  modified  power  law  (MFL)  in  reduced  time.  The  coefficients  of  the  MPL  are 
found  through  nonlinear  inversion  of  the  data  using  a  Marquardt-Levenberg  algorithm.  Time-temperature  superposition  is 
employed  and  the  horizontal  shift  function  for  each  strain  level  is  given  by  an  Arrhenius  expression.  The  propellant  is 
nonlinearly  viscoelastic  since  the  logarithm  of  the  vertical  ^ift  function  is  a  quadratic  function  of  the  logarithmic  strain. 
Theoretical  predictions  of  time-dependent  stress  versus  time,  failure  stress  versus  failure  time  artd  failure  stress  versus  strain 
rate  quantitiatively  agree  with  the  experimental  behavior  of  the  propellant  determined  at  various  strain  rates  and 
temperatures. 
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1.  INTRODUCTION 


The  occurrrace  of  high-amplitude  i»essuie  waves  in  guns  may  lead  to  catastrc^rtik;  breech 
blows.  Several  factms  duit  cause  the  inessure  waves  include  a  malfunctioning  igniter,  reduced  charge 
permeability,  a  high  initial  gas  generation  rate,  and  charge  compaction  which  results  in  grain  frachire 
(MayandHorst  1979;  Keller  and  Horst  1989).  Ckain  fracture  leads  to  an  increase  in  the  mass  burning 
rate  of  the  propelling  charge  through  an  increase  in  the  surface  area  available  fm*  combustion.  The 
growth  and  coalescence  of  stress-induced  microfractures  alter  the  microstructure  of  gun  prc^Uant  and 
can  dramatically  increase  the  apparent  burning  rate  in  damaged  propellant  relative  to  undamaged 
propellant  (Gazonas  et  al.  1991). 

The  dynamics  of  crack  growth  in  granular  propellant  are  critically  dependent  upon  the  crack 
velocity  relative  to  the  propellant  bum  rate.  Slowly  growing  microcracks  that  propagate  as  a  result  of 
intergranular  stress  in  the  propellant  bed  may  be  entirely  consumed  by  rapid  burning  of  the  propellant 
grains.  In  some  very  high  bum  rate  (VHBR)  hivelite-based  propellants,  proposed  for  use  in  traveling 
charges,  the  rq)patent  burning  rate  approaches  300  m/sec  (May  et  al.  1984);  a  surface-fracturing 
mechanism  was  proposed  to  explain  the  npid  burning  rates  observed  in  VHBR  propellants.  The  burning 
rates  of  conventional  gun  propellants  are  considerably  less  than  in  VHBR  propellants.  For  example,  the 
burning  rates  of  single  grains  of  undamaged  M30  propellant  are  pressure-dependent  and  approach  10 
cm/secat60MPahydro^aticpressure(Gazonasetal.  1991).  Alternatively,  rapidly  growing  microcracks 
thatare  hydraulically  driven  by  gaseouscombustion  products  (Wuetal.  1992;  Griffiths  and  Nilson  1992) 
might  accelerate  the  ipparent  burning  rate  through  convective  heating  of  the  propellant  Cracks  in  many 
materials  propagate  at  velocities  which  range  from  quasistatic  (Atkinson  and  Meredith  1 987)  to  dynamic 
(Swanson  1982);  the  upper  bound  on  crack  velocity  is  theoretically  limited  by  the  Rayleigh  wave 
velocity,  under  usual  types  of  loading  in  isotropic  elastic  media  (Sih  1968).  However,  supersonic  crack 
velocities  have  been  observed  in  solids  subjected  to  strong  laser  pulses  (Winkler  1973)  or  high  intensity 
electron  beams  (Cherepanov  and  Borzykh  1993).  Sbearwave^ieedsinlow-densitypropellantmaterials 
are  approximately  2  km/sec  (Costantino  1983),  which  indicates  that  maximum  crack  velocities  in 
conventional  propellants  are  about  four  orders  of  magnitude  greater  than  propellant  burning  rates. 


Hence,  for  most  interior  ballistic  applications  that  involve  high-amplitude  fuessure  wave  loading  of  the 
propellant,  microcrack  growth  rates  will  far  exceed  the  burning  rate  of  the  propellant  The  di^uuity  in 
tinic  scales  between  physical  processes  (crack  speed  and  burning  rate)  is  further  exacerbated  by  the  fact 
that  not  ail  regions  of  the  granular  charge  are  ignited  simultaneously;  the  flamespreading  rate  depends 
upon  the  gas  permeability  of  the  granular  propellant  charge.  A  consideration  of  these  observations  could 
simplify  interior  ballistic  numerical  model  development  because  in  some  instances,  the  mechanics 
problem  can  be  temporally  decoupled  from  the  combustion  problem. 

However,  the  presence  of  microfractuies  in  gun  propellant  is  not  necessarily  deleterious  to  the 
mechanical,  combustion,  or  ~  Inerability  response  of  the  propellant;  it  is  known,  for  example,  that 
distributed  microcracks  in  certain  composite  materials  can  improve  their  toughness  ot  global  resistance 
to  fracture.  In  principle  then,  a  gun  propellant  could  be  designed  to  "programmably"  fracture  in  a 
predictable  and  reproducible  manner.  This  proposal  is  contrary  to  traditional  strategies  of  propellant 
formulation  that  aim  to  minimize  the  fracturability  of  the  propellant  through  the  addition  of  energetic 
plasticizers  and  other  ingredients.  The  microstructural  fabric  of  the  propellant  also  plays  a  major  role 
in  subsequent  growth  and  coalescence  of  microcracks.  For  example,  Fong  and  Moy  (1982)  suggest  that 
propellants  that  contain  ground  spherical  nitroguanidine  oxidizer  will  have  fewer  possible  crack 
propagation  sites  than  propellants  that  contain  long  needle-like  nitroguanidine  crystallites.  Grain  size, 
orientation,  pore  size,  shape  and  distribution,  and  a  multitude  of  other  structural  elements  that  compose 
the  microstructural  fabric  of  the  propellant  will  undoubtedly  influence  the  initial  distribution  and 
subseqirent  growth  and  coalescence  of  microcracks.  The  proposed  scheme  recognizes  that  microcracks 
are  present  in  nearly  all  manufactured  materials  and  that  the  effects  of  microcrack  growth  can  be 
predicted  and  therefore  advantageously  used  in  the  design  of  new  propellants. 

The  idea  that  gun  performance  is  improved  if  additional  surface  area  is  introduced  into  the 
propellant  charge  is  not  new.  For  quite  some  time,  ammunition  designers  have  used  multi-perforation 
grain  design  concepts  to  alter  the  mass  generation  rate  of  propellant  during  burning  (for  example,  see 
Military  Exploaves  ( 19SS),  and  Engii^ring  Design  Handbook  ( 1 965)).  Multiperforation  grain  designs 
allow  for  the  "progressive"  increase  in  the  total  surface  area  of  the  propellant  grain  as  it  bums.  By 
mathematically  modeling  the  progressive  growth  of  stress-induced  microfractures,  one  could,  in 
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principle,  design  high-performance  yet  "safe"  propellant  charges  far  fracture. 


An  essential  feature  of  the  strategy  just  described  is  the  development  of  a  mathematical  model 
that  can  predict  the  incremental  time-dependent  increase  in  stress-induced  fracture  surface  area  in  the 
propellant  at  temperature  and  loading  rate  conditions  relevant  to  the  gun  environment  Recent  advances 
in  continuum  mechanics  make  it  possible  to  characterize  and  quantify  microcrack  growth  using  concepts 
from  the  relatively  new  field  of  continuum  damage  mechanics  (CDM).  Krajcinovic  (1989)  provides  an 
^n-depth  review  of  CDM  from  both  microphysical  and  phenomenological  standpoints.  The  CDM  model 
predictions  of  fracture  surface  area  could  then  be  quantitatively  compared  to  the  progressive  develop¬ 
ment  of  fracture  surface  area  in  the  propellant  derived  by  conducting  closed  bomb  tests  on  incrementally 
damaged  propellant.  Progress  with  this  concept  began  in  a  study  that  examined  the  effect  of  various 
mechanical  variables  such  as  strain,  strain  rate,  and  temperature  on  the  combustion  response  of  two 
conventional  gun  propellants  (Gazonasetal.  1991).  In  that  study,  it  was  found  thatsurface  areaevolution 
in  M30  and  JA2  gun  propellants  is  relatively  insen^tive  to  strain  rate  when  compared  to  the  effects  of 
temperature  and  the  degree  of  axial  deformation. 

In  this  report,  a  uniaxial  nonlinear  thermoviscoelastic  constitutive  model  with  time-dependent 
damage  is  developed  for  the  conventional  tank  gun  propellant,  M30.  This  work  extends  previous  woik 
that  modeled  loading  rate  behavior  of  the  propellant  over  nearly  five  decades  in  strain  rate  (Gazonas 
1992, 1993).  Wave  speeds  in  porous  propellant  media  have  been  directly  measured  as  a  function  of 
hydrostatic  pressure  (Costantino  1983)  and  estimated  from  propellant  bed  compaction  studies  (Kooker 
1990).  However,  constitutive  models  for  gun  propellants  are  virtually  nonexistent  This  work  represents 
a  continuing  effort  to  mathematically  model  the  constitutive  response  of  solid  propellant  necessary  in 
two-phase  (solid  and  gas)  flow  numerical  models  of  interior  ballistic  processes.  Furthermore,  the  work 
offers  a  means  for  mathematically  modeling  the  progressive  growth  of  microcracks  that  could  be  used 
in  designing  new  propellants  with  the  expanded  "firacturable"  propellant  design  concept  described 
previously. 

This  report  employs  a  uniaxial  specialization  of  a  general  three-dimensional  constitutive  theory 
for  nonlinear  viscoelastic  materials  with  damage  (Schapery  1981).  Herein,  the  previous  rate-dependent 
model  is  extended  to  incorporate  temperature  effects  by  replacing  real  time  variables  with  reduced  time 


3 


aiid  frilnie  alien  vmnsatniii  file  qoaiititaiivefy  agree  experime^ 

deinmlfied  at  vtfion  str^  ratea  and  tempentitfes. 


2.  THE  CONSTITUTIVE  THECMIY 

Biot  (1954)  developed  the  unified  ttieory  of  linear  viscoelasticity  based  upcm  the  thewy  of 
irreversible  thermodynamics.  Later,  the  themy  was  extended  to  include  nonlinear  thermoviscoelastic 
media  (Schi^)^  1964)  md  distributed  damage  effects  using  a  so-called  "modified  superposition 
imegral"  (Schiqtery  1981).  The  uniaxial  strain  e  in  a  material  subjected  to  a  uniaxial  stress  a  can  be 
written  with  the  notation  of  Schapery  (1989)  as. 


t 

e(t)  *  Ej  J  D(t-T)jlf  dr  .  (1) 

0*  dt 

The  integral  in  (1 )  is  also  known  as  an  hereditaiy  or  ccmvolution  of  die  functions  D  and  f.  The  strain  is 
thus  seen  to  be  dependent  uptm  the  oitire  hiaoryttf  appUcadon  (rf^stness  and  not  just  (m  the  instantaneous 
value  as  in  the  case  of  nonlineaily  elastic  matmials.  It  is  assumed  that  e  s  o  s  f  s  0whent<0and 
D(t-x)  s  0  when  t  <  t.  Furthermore,  the  lower  limit  in  die  integral  is  O'  to  allow  for  a  possible  jump 
discontinuity  in  the  function  f(t)  at  t  s  0.  All  subsequent  hereditaiy  integrals  iqqiearing  with  the  lower 
UmitatOaie  to  be  interpreted  as  O'  as  well  D  is  the  uniaxial  creep  compliance,  E^  is  aconstant  reference 
modulus  with  units  of  stress,  andf  is  given  as  a  product  of  two  functions  g,  and  g,,  as. 


f=f(o.SJ*g,«j)  g,„(S^  ,  (2) 


in  which  g,  is  typically  written  as  a  power-law  funcdon  of  stress  and  g,^  is  a  funcdon  of  the  damage 
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parameter.  The  are  k  damage  parameters  used  to  describe  microstruchiral  changes  that  occur  which 
alter  the  material's  fabric.  In  previous  work  (Gazonas  1992.1993)  it  was  shown  that  a  single  scalar 
damage  parameter,  S^,  was  sufficient  for  modeling  mkrocrack  growth  in  M30  gun  propellant.  One  could 
argue  that  a  one-dimensional  model  cannot  further  our  understanding  of  thiee-dimensioiial  microi^ysi- 
cal  fracture  processes  that  occur  in  a  propellant  grain.  Nevertheless,  the  simplification  i4>pears  to  be  valid 
for  describing  the  global  constitutive  behavior  of  M30  propellant  and  a  variety  of  other  materials 
including  marine  sediment,  composites,  polymers,  rocket  propellant,  and  ice.  Regardless  of  the 
dimensionality  of  the  proposed  model,  it  should  ultimately  be  reducible  to  a  one-dimensional  form  diat 
can  be  verified  through  experimentation.  If  a  sufficient  number  of  material  tests  are  conducted  using 
deformation  histories  similar  to  those  experienced  by  the  gun  prqiellant  during  Hring,  one  can  develop 
empirical  models  thus  circumventing  the  need  for  the  development  of  more  complex  microphysical 
4nodels. 

The  damage  parameter  in  (2)  is  stress-history  dependent  and  is  expressed  as. 


in  which  a,  and  q  are  positive  constants,  and  1 1  denotes  the  absolute  value  of  the  quantity.  The  "crack 
tip  material-related  coefficient"  f,  can  depend  on  time  and  temperature  and  material  aging  effects 
(Schapery  1981,1989).  For  example,  extension  of  the  isothermal  theory  to  account  for  transient 
temperature  effects  involves  replacing  f,  in  (3)  by  the  factor  (1/a,.)  and  using  reduced  time  in  Eqns.  (1) 
through  (3).  The  thermomechanical  shift  factor,  a,. ,  is  related  to  entropy  production  in  Schapery's 
thermodynamic  theory.  The  extension  of  the  model  to  incorporate  thermoviscoelastic  behavior  of  the 
propellant  is  provided  in  a  subsequent  section. 

The  underlying  physical  basis  for  the  constitutive  theory  with  distributed  damage  is  a  crack 
growth  "law"  that  relates  the  speed  a  of  a  single  crack  to  a  power-law  function  of  the  mode  I  stress 
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iiurasity  factor,  a  =  A  K,.  in  which  A«idQare  empirically  d^nnined  constants.  SchiiMfy  (1975) 

arrives  at  this  result  by  generalizing  Barenbliut's  (1962)  cohesive  ciack  tip  model  for  linear  elastic  media 
to  linearly  viscoelastic  media.  The  crack  propagation  theory  is  later  generalized  (Schapery  1981,1984) 
to  encompass  nonlinear  viscoelastic  behavior  through  the  use  of  the  so^alled  "pseudo- variables.”  The 
use  of  pseudo-variables  for  describing  damage  effects  in  nonlinearly  viscoelastic  media  effectively 
reduces  the  problem  to  one  of  nonlinear  elasticity  since  stiess-pseudostrain  plots  are  single-valued 
(Schapery  1982).  In  an  analogous  fashion,  a  pseudostiess  transformation  could  be  invoked  to  create  a 
single-valued  curve  in  woiksoftening  materials. 

A  complete  derivation  that  extends  the  single-crack  theory  to  describe  the  global  response  of 
the  medium  due  to  distributed  microcracking  is  beyond  the  scope  of  this  report  but  is  addressed  in 
Schapery  (1981)  (Equations  122  through  126).  The  derivation  uses  a  relation  between  crack  speed  and 
a  power-law  function  of  a  generalized  J-integral  (Schapery  1981, 1984). 

2. 1  Strain-History  Inputs.  Equations  (1)  through  (3)  are  suitable  for  characterizing  a  material's 
strain  response  if  stress  is  a  controlled  input  fcMr  the  test  However,  if  strain  is  a  controlled  input  for  the 
test,  then  these  equations  can  be  inverted  to  predict  stress  as  a  fiuiction  of  strain  history, 

o(t)  =  Oj  e“'*  g^(Sp  sgn(€?)  .  (4) 

in  which  sgn  is  the  signum  function,  o,  is  a  constant  parameter  with  units  of  stress, 

'(Schapery  1989),  and  the  pseudo-strain,  e°,  is  related  to  strain  history  by. 


inwhkdiECOistlieielaxationiDodidiifc  AdaiiMifep«nm^S,for8tndBU«acyiB|iiitt(Sciiq)e^  1989) 
is 


S. 


fidt  . 


0 


(6) 


in  which  o,,  Oj.  q.  and  r  are  mateiial-depeixlent  cinstants. 

2.2  Extension  to  ThermoviacoelMticitv  with  Dmaye  TheexleiiaiottoflteiBoddtoiiicorp^ 
thermoviscoelasticity  widi  damage  is  pofonned  by  replacing  die  tune  variable  t  widi  a  reduced  time 
variable^  in  the  previous  equations.  Biiriieimme.detennination  of  die  mMlerrelaxatioa  modulus  E(4) 
is  tiadidonally  obtained  by  perfrmning  stress-relaxation  tests  ri  various  tenqientines  and  then  rigidly 
shiftily  the  time  scales  fmr  each  of  die  tests  to  f<xm  a  continuous,  single- valued  carve  (Morland  and  Lee 
1960).  The  single-valued  curve  dien  refuesents  the  relaxadmi  modulus,  (or  creep  compliance  if  a  unit 
Ifeavisidestiess-hisUMy  is  iqiplied)  at  the  reference  temperature,  for  very  short,  intermediate,  and  very 
kMigdmes.  If  the  master  relaxadrui  modulus  can  be  fcamed  through  horizontal  translations  of  the  data, 
the  rrurterial  is  termed  thermmheologically  ample.  The  master  relaxadmi  modulus  is  plotted  as  a 
function  of  a  reduced  time  parametn' 4  defined  by. 


(7) 


in  which  a,,  is  the  thermomechanical  shift  function  defined  as  a,  si  at  the  rtference  temperature.  The 
prediction  of  time-dependent  stress  in  (4)  is  then  made  by  replacing  real  time  in  (5)  with  reduced  time, 
which  yields  for  die  pseudostrain. 
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8ectk)ii.  it  is  shown  how  coostiitt  flnifl  itte  coi^vetfioii  tests,  cowhicied  at  vaikMis 
be  used  todeteimine  E(4)- 


3.  EXPERIh«ENTAL  RESULTS 

Experimental  results  fw  isothennsi  unisxial  comfKessioa  of  M30  pn^Unnt  St  two  stndn  files, 
0.01  and  210  sec‘‘,  and  three  temperatures,  -30*.  21* ,  and  60*  C  ^ipear  in  Figure  1.  Eadi  siress-tiroe 
curve  is  a  composite  curve  formed  from  die  average  oi  five  tests.  Tlie  uniaxial  compression  tests  are 
performed  using  an  MTS  hi^  rate  8 1 0  material  lest  system  described  in  more  detail  elsewhere  (Gazonas 
1991).  The  maximum  piston  velocity  is  approximately  12  m/sec,  wUdi  limits  die  axial  compressive 
strain  rate  to 500 sec‘(Mi2S-mmq)ecimens.  Constantstrainraletestsareperformedbycompulercontrol 
of  the  {ristcMi  velocity  dirough  feedback  fiom  an  externally  mounted  linear-variable  differential 
transformer  (LVDT,  MTS  Model  244.1 1).  FOice  is  measured  with  a  60-kN  quartz-piezoelectric  force 
transducer  (Kisder  Type  9031 A)  that  is  mountnl  cm  die  moving  pisttm  (Figure  2).  Specimen 

di^lacements  ate  ccniectedfcv  apparatus  di^ortkmdiat  has  a  ineasured  stiffness  of  abom  92  kN/mm. 
Specimen  stiffness  ranges  from  4.5  to  14  kN/mm  id  room  temperature  fcv  the  strain  rate  range 
investigated,  ^lecimens  are  diermally  conditioned  at  temperature  for  about  1  hour  before  testing  in  a 
Thermotron  conditicming  oven/refrigerator  (Controller  Model  5200)  that  houses  both  iqqier  and  lower 
pistons.  The  chemical  composition  and  nominal  specimen  dimensions  of  M30  aiqiear  in  Table  1. 


4.  MASTER  RELAXATION  MODULUS  DETERMINATION 

The  master  relaxation  modulus  is  experimentally  determined  using  the  results  of  constant  strain 
rale  tests  conducted  at -30*,  21*,  and  60*  C.  In  eartierworic,  it  was  shown  diat  die  relaxation  modulus 
f(MrM30,JA2,aiidXM39propellantsisderivaUefromconstant8tninrmetestdata(Qazonas  1991).  For 
a  ccmstant  strain  rate  test,  the  one-dimensional  linear  viscoelastic  ccmstitutive  equatkm  widKMit  damage 
is  written  as. 
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t 

0(0  •  f  dt  0) 

0  dc 


i^acoostaotstnianaeiiqpiithistary.le.,  e«  e  t.  (9)becoines, 


t 

<^0  >  e  J  E(t-t)  dr  .  (10) 

0 


The  secant  fflodohis  is  defined  as. 


E  «  A 
•  e  • 


(11) 


and  since  e  «  e  t. 


E  >  t 


-l 


% 

J  E(tt)do 


(12) 


Taking  doivatives  of  bodi  sides  (12)  lesolts  man  expression  for  the  rdaxalioo  moduhis. 


dE 

B(t)  »  ^  +  t 


(13) 


Eqmdion  (13)  can  be  rewrittra  as 


E(0-E,(1. 


(14) 
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IMi«  •  inpoAgr  of  loprittM  aA.  ii  10  4(Iq|^  «  (14)  CM  te  lewri^ 


m  -  B,(i.  ♦ 


^fiA 

"3517 


(15) 


If  the  slope  (rf  the  lpg|^,  versus  log^^t  curve  is  mudi  less  dun  unhy,  then  die  rdaxadoo  modulns  is 
essentially  e^valentto  the  secsntmodulus.  TodeiermineiftheielaxationinodulosofMSOisnoiilhiear 
(strain  dependent),  one  can  plot  the  secant  modulus  ate  at  stndn  levels  of  2. 5, 10. 15.20,25,30.  and 
35%  versos  lethioed  time  4  (Rsme  3). 


If  the  thermomechanical  shift  fimcdon  a,  dqiends  primarily  on  ten^pnature  and  not  on  sindn 
levelin(7),thenaf  aaf(D.  Purtheimote.  if  tea^ieratuie  is  timewise  constant,  Equadon  (7)  reduces  to 
^st^a^  the  effects  <tftrandenttempentoiedianges  in  the  gun  environment  may  be  negligiMe  relative 
to  effects  (tf  the  propagation  of  mechanical  disturbances  (P  and  S  waves).  The  therriiomedianical  shift 
function  is  described  by  an  Arrhenius  expression. 


(16) 


in  ududiQis  die  acttvationenetgyODcal/taKde),RisdMimivet8algasconstant  (1.987  cal/g>nKdeK),  and 
T,  is  the  reference  temperature  (K).  The  activation  enei^  is  27.4  kcai/ioBote  determined  from  the  slope 
of  die  log, 0  (Sf)  versus  1/T  curve  since  d(lQg,Q(a,))/(KlfT) « (y2J03R  (Hgure  4).  Aldioogh  die  data  are 
Umiied,  it  appears  dud  the  dremromechanical  ddft  fdncdon  exhibils  some  strain  dqiendmioe.  The 
amount  of  strain-dqpendent  shift  is  greater  at -30*  C  than  at  60*C(Rgure  4).  From  die  stress  relaxation 
dmattfUeband  Leadme  (1992),  an  acdvadon  oiersy  of  29.6kcayhiole  is  derived  that  is  only  8%  greater 
than  die  activation  energy  repmted  hoein;  die  difference  in  the  calculated  activation  energy  may  also 
be  attributed  to  the  fiKt  dud  the  stress  relaxadoo  experiments  Lieb  and  Leadore  ( 1992)  are  conducted 
on  a  different  kd  of  mnldpetfoiaied  M30  propellant,  whoeas  die  test  remits  repotted  herein  are 
conducted  on  atdid  stidc  M30  propeOard. 
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Since  die  material  is  nonUneariy  viscoelastic,  a  vertical  shift  function,  a^,  is  used  to  colla(we  the 
relaxation  moduli  determined  at  each  strain  level  (H^ure  3)  onto  a  single  master  curve.  The  logarithm 
of  the  vertical  shift  function  is  empirically  determined  to  beaquadiatic  fiinctiMi  of  the  logarithm  of  strain 
(Figure  5).  The  amount  of  vertical  shift  is  unity  when  the  strain  level  e  equals  die  reference  strain  level 
e,  of  2%.  Figured  illustratesan  overiaycrthe  data  and  the  master  relaxatkm  modulus  (at2%  strain)  which 
takes  the  fmm  of  a  modified  power-law  (MPL)  in  reduced  time;  the  coefiicients  of  the  ate 
determined  using  a  nonlinear  least  squares  inversion  technkjue  (Marquardt-Levenberg  algorithm).  The 
MPL  is  written  as. 


(Eq-eJ 

(1.  +  f )“ 


(17) 


in  which  »  E(0^)  and  E. = E(««)  ate  the  shrut-term  and  long-term  asymptotes  of  the  relaxation  modulus, 
and  t,  and  n  are  additional  firee  constants.  The  constants  assume  values,  E,  s  4,470  ±  628  MPa,  E^  %  453 
±94MPa.n-0.171  ±.031,andt,=:4.6x  10^±5.2x  lO^secs.  Equation  (17)  increases  the  range  of 
behavior  of  M30  over  a  broad  time  plectrum.  For  example,  if  the  propellant  is  strained  to  2%  at  room 
temperature,  the  secant  modulus  relaxes  from  4,470  MPa  to  453  MPa  (about  an  order  of  magnitude)  in 
about  28  hours.  The  same  relaxation  requires  about  1 .6  minutes  at  60°  C  and  about  6.4  years  at  -30°  C. 

Substitution  of  the  relaxation  modulus  E(4)  into  (8)  and  then  substitution  of  the  pseudostrain  into 
(6)  and  (4)  completes  the  extenrion  of  the  model  to  thermoviscoelasticity.  If  the  theory  is  cturect,  the 

I  2 

form  of  the  damage  function  gjgfSg)  =  lO'*’’  -  T  i  -  ^  identical  to  that  determined  previously  for 

M30  (Gazonas  1991,1992),  in  which  %  =  logjo(S^)/log,g(S^).  The  assumption  is  that  all  temperature 
effects  are  introduced  into  the  damage  function  g^fSp  and  damage  parameter  S,  through  the  |»eudostrain 
in  (8)  which  is  a  function  of  temperature-reduced  time. 

How  well  the  model  fits  the  dataforapaiticular  material  will  depend  upon  whetherasingle  scalar 
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PievkMtwoA  accontdy  oiodd^kMtfiig  nie  dfecift  (over  oearty  fiwe  decades  ill  stnia  rale)  in  M30 
with  a  single,  scalar  damafe  parameter  (Hgure  7).  bi  geneial.  multipk  damage  paraBMsters  ate  specified 
using  (2);  a  damage-reduced  time  duft  function  approach  might  also  be  employed  to  accoiutt  for  mcue 
complex  temperature-dependent  behavior. 

The  currentmodel  is  furdiersimplifiedifapower-law  relaxation  modulus(Gazonas  1992, 1993), 

Edl^Eit-  .  (18) 

is  used  rather  than  the  MI^  relaxation  modulus  in  (17).  Replacing  real  time  with  reduced  time  in  (18) 
and  then  substituting  (18)  into  (8)  with  %  =  t/a^  E,  =  E^  and  for  a  constant  strain  rate  input  histtny. 

e  s  e  t ,  yields  an  expression  for  the  pseudostrain. 


.»-a 


t^t)  *  Sf  e  I 


(19) 


Substitution  of  (19)  into  (6)  with  f^  =  (l/a^)  and  Oj  =  yields  for  the  damage  parameter. 


t  1 

•r  (i-n>l  +  l 


(20) 


Time-dependent  stress  (Figure  7)  is  predicted  by  substituting  (19)  and  (20)  into  (4)  using  the  ellipsoidal 
damage  function  g„(S^  and  constants  detnmined  previously  ((3azonas  1992, 1993).  Model  predictions 
of  time-dependent  stress  agree  faiily  well  with  the  experimental  data  (Figure  8). 

The  failure  stress  (maximum  stress)  in  M30  propellant  decreases  as  both  failure  time  (time 
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required  to  reach  maximum  stress)  and  temperature  increases.  The  coefficient  of  die  power  law  relation 
is  about  •  0.10  for  tests  conducted  at  -30”,  21”.  and  60”  C  (Hgure  9).  The  failure  stress  at  a  given 
temperature  increases  by  a  factor  ranging  from  2  to  3  over  nearly  five  decades  in  ^rain  r^.  Furthermore, 
the  mechanical  behavior  of  M30  is  dominated  by  temperature  changes  since  the  failure  stress  at  a  given 
strain  rate  also  increases  by  a  factor  ranging  from  2  to  3  when  the  temperature  increases  fircMn  21”  C  to 
60”C;  equivalent  changes  in  the  magnitude  of  the  failure  stress  in  the  {vq[)ellant  are  seen  by  changing 
strain  rate  by  nearly  five  orders  of  magnitude.  Failure  strains  are  seen  to  increase  slightly  with 
temperature:  however  at  a  given  temperature,  failure  strains  are  insensitive  to  strain  rate.  This  is  seen 
by  substituting  the  power-law  expression  for  failure  time  and  strain  rate  (Figures  9  and  10)  into  the 

expression,  =  t,  e,  and  recognizing  that  the  result  is  a  constant 


5.  SUMMARY  AND  CONCLUSIONS 

A  nonlinear  theory  of  viscoelasticity  with  damage  (Schtqiery  1 98 1 )  has  been  extended  to  account 
for  the  effects  of  both  temperature  (-30”  m  60”  C)  and  strain  rate  ( 1 sec'  to 400  sec')  on  the  constitutive 
response  of  M30  gun  propellant  Theoretical  time-dependent  stresses  in  the  propellant  subjected  to  a 
timewise  constant  temperature  field  and  arbitrary  input  strain-history  are  obtained  by  replacing  real  time 
variables  with  reduced  time  variables.  The  theoretical  predictions  of  time-dependent  stress  versus  time, 
failure  stress  versus  failure  time,  and  failure  stress  versus  strain  rate  quantitatively  agree  with  the 
experimental  behavior  of  the  propellant  determined  at  various  strain  rates  and  temperatures. 

The  thermomechanical  shift  function  for  M30  propellant  is  determined  by  shifting  data  derived 
from  uniaxial  constant  strain  rate  compression  tests.  The  shift  function  follows  an  Arrhenius  law  with 
an  activation  energy  of  approximately  27.4  kcal/mole.  This  value  compares  well  with  an  activation 
energy  of  29.6  kcal/mole  derived  directly  from  the  stress  relaxation  test  results  of  Lieb  and  Leadore 
(1992). 


The  relaxation  modulus  of  M30  propellant  is  nonlinear  (strain  dependent),  and  the  logarithm  of 
the  vertical  shift  function  is  a  quadratic  function  of  the  lo^irithmic  strain.  The  relaxation  modulus  for 
M30  propellant  has  been  obtained  in  the  form  of  a  modified  power-law  in  reduced  time;  the  coefficients 
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of  die  IfffL  aie  fooMi  dmogb  aotfifiear  tevcfskn  (rf  tte  dWa  oiiaf  t  Maiqiiwdt^jev^^ 

However,  die  aiudyeb  become  sofloeiidiatsta^^  power-lew  foim  for  die  fdaxatkmfiiodidiis 

detennifiedinearti«rwaric(Oazoiiasl992. 1993)aiidiiicorpofitiiigalliionliiiearitiesaiiddem«feindie 
function  f(a.S^) )  (Equation  2). 

The  effects  of  tempoatuiedcmiinate  strain  rate  effects  in  M30|m>pellant  as  the  magnitude  of  the 
failure  stress  (maximum  stress)  is  equally  affected  by  achange  in  temperature  from  21°  to  60”  C  ((M’21° 
to-30”  C)  or  by  a  change  in  strain  rate  that  ^lans  nearly  five  (Ntlersttf  magnitude.  The  master  curve  is 
also  significant  nnce  measurements  made  at  labtvatmy  experiment  time  scales  are  extrapolatable  to  very 
sh(»t  or  very  long  times.  For  example,  if  die  |»opellant  is  strained  to  2%  at  romn  temperature,  the  secant 
modulus  relaxes  from  4,470  MPa  to  4S3  MPa  (about  an  order  of  magnitude)  in  about  28  hours.  The  same 
relaxation  requires  about  1.6  minutes  at  60”  C  and  about  6.4  years  at  -30”  C.  Most  interim  ballistic 
applicaticHis  require  extrapolation  of  material  (Koperties  to  very  short  times,  on  the  cmler  of  milliseconds. 
The  model  developed  herein,  however,  is  also  useful  for  extrapolating  material  properties  and  constitu¬ 
tive  behavior  to  very  long  times  and  could  be  useful  in  the  long-term  design  of  storage  facilities  for 
propellants.  Other  factors  such  as  the  effects  of  humidity  and  aging  effects  of  gun  propellants  are  also 
important,  but  not  well  understood,  and  require  further  testing  to  aid  in  model  development  and 
validation. 
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Table  1.  Chemical  Compositkn  and  Nomiaid  Specimen  Dimensioiis 

(rf  M30  Gun  Propellam 


Gnnpoiient 

% 

Nitrocellulose 

28.0 

12.6  %NC  Nitration 

Nitroglycerin 

22.0 

Nitroguanidine 

48.0 

Ethyl  Centralite 

2.0 

100.0 

Length  (mm) 

25.4 

Diameter  (mm) 

12.3 

IS 
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Figure  7.  Overiav  of  Observed  (symbols^  and  Theoretical  (solid  lines)  Stress  versug  Tune 
Curyes  for  M30  Gim  Propallaat  Under  Constant  Strain  Rate  Deformation. 
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Figure  8.  Overlay  of  Observed  (symbols)  and  Theoretical  (8oUd  lines)  Stress  versus  Time 
Curves  for  M30  Gun  Propellant  at  Two  Strain  Rates  fO.Ol  and  210  ae<r»l  and 
Temneratures  f.30>  and  (Sft»  V 
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AMSRL-Wr-WA, 

R  Rogers  (394) 

B.  Moore  (394) 

A.  Boan  (394) 


AMSRL>WT-WB. 

F.  Brandon  (120) 
W.  D'Amico  (120) 
AMSRL-WT-WC. 

J.  Rocchio  (120) 
AMSRL-WT-WD. 


A.  ZieUnski  (120) 
J.  FOweU  (120) 
AMSRL-WT-WE, 


J.  TemfMiley  (120) 
J.  Thomas  (394) 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

ThU  Labompiy  undeittkei  «  continMii^  cflbn  to  improve  the  quality  of  the  repoitt  it  pribiiRiet.  Your 
oommenii^niwen  lo  die  UemMliiefltkMif  bdow  will  aid  us  in  our  eflbfts. 

1.  ARLReooit Number  AHr-TK-469 _ DmerfRepoit  9A _ 

2.  Dale  Report  Received _ 

3.  Does  dds  report  sadafy  a  need?  (Cnmment  oo  puipoae,  idased  pioject,  or  odier  area  of  imerest  fcr 

wldcli  die  report  win  be  used.) _ 


4.  SpedficaUy,  how  is  die  report  being  used?  (Infimnadon  soune.  design  data,  procedure,  source  of 
ideas,  etc.) _ 


S.  Has  the  information  in  this  report  led  to  any  quarttitadve  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  fdease  elaborate. _ 


6.  General  Cmnments.  What  do  you  think  should  be  dianged  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ _ 

Street  or  P.O.  Box  No. 


Qty,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  (riease  provide  the  Current  or  Correct  address 
above  and  die  CBd  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  shed,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


Department  or  the  Army 


OFROAL  BUSINESS 


BUSDOESSREPLYMAIL 

BKraiBsiairiiinii.iiFt,iB 

Pesiagc  mitt  be  aanS  by  addrtuw. 


Director 

U^.  Army  Research  Laboratory 
ATTN:  AMSRLOP-a-B  (Tech  Ub) 
'Aberdeen  Proving  Ground.  MD  21005-5066 


